The discovery of new warm debris disks around F-type stars. 
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ABSTRACT 



We report the discovery of four rare debris disks with warm excesses around F-stars, 
significantly increasing the number of such systems known in the solar neighborhood. 
Three of the disks are consistent with the predictions of steady state planetesimal disk 
evolution models. The oldest source, HD 169666, displays a dust fractional luminosity 
too high to be in steady state and we suggest that this system recently underwent a 
transient event of dust production. In addition, two spectra of this star separated by 
~three years show silicate emission features, indicative of sub micron- to micron-sized 
grains. We argue that such small grains would be rapidly depleted and their presence 
in both spectra suggests that the production of small dust is continuous over at least 
on few years timescale. We predict that systems showing variable mid-infrared spectra. 
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if exist, will provide valuable help in distinguishing the possible scenarios proposed for 
dust replenishment. 

Subject headings: circumstellar matter — infrared: stars 



1. Introduction 



Many main-sequence stars exhibit excess emission at mid-to-far infrared (IR) wavelengths, 
implying the existence of a ring or disk of circumstellar dust grains, which re-radiate the absorbed 
starlight at longer wavelengths. Due to mutual collisions, radiation pressure and the Poynting- 
Robertson (PR) drag, the lifetime of the emitting dust grains is much shorter than the age of the 
star, and the dust particles are thought to be contin uously replenished by the erosion of larger 



planetesimals ([Backman Sz Paresce 



1993 



WyattI I2OO8I ) . Thus, the existence of a dust 'debris' disk 



implies the presence of a planetesimal belt around the star. 

In most known debris disks around F, G, and K-type stars the dust grains are relatively 
cold (50-90 K), s i milar in temperature to the dust located in the Kuiper Belt in our solar system 
(jCarpenter et al.ll2009l ). Interestingly, warm debris disks - where the emission (in F^) peaks at 
A > 40 /im, corresponding to T^ust ^ 125 K - are rare. If the infrared emission in these warm 
disks originated from blackbody grains, they would be located at a radius of 1-10 AU. Thus, in 
'solar-system terminology' these disks are either similar to the main asteroid belt or are situated 
between the main belt and the Kuiper Belt. 

So far ~20 warm debris disks have been discovered around FGK-type stars. Most of them en- 
circle young stars, that are often me mbers of known young open clusters or moving groups typically 



located at larger ( > 100 pc) distances (IChen et al 



Hines et al. 



2006 



Lisse et al. 



2008 



Rhee et al 



2005 



2008 



Smith et al 



Currie et al.ll2007l . I2OO8I : ICxOrlova et al.ll2007l : 



20081 ). In these young (<100 Myr) 



systems the source of warm debris dust in the region of terrestrial planets may be an outwards 



propagating ring of 



teenyon &: Bromley 



jlanetesimal formation and /or the products of massive protoplanet collisions 



2004l : lMever et al 



around older solar-type stars (jSong et al..l 



20081) . A smaller number of warm disks have been identified 



2005 



Beichman et al 



20051 . 



2006 



Wvatt et al.ll2005l l. Ah 



these old systems exhibit high fractional luminosity (7;^ > 10~^), that cannot be explained with a 
steady-state as teroid belt evolutio n model, where the planetesimals are co-located with the warm 
dust. Instead, IWvatt et al.l (|2007l ) suggested that these systems are in a transient state. Recent 
collisions between large asteroids or the erosion/sublimation of planetesimals scattered from an 
outer reservoir into the in ner regions due to a dynamical insta bility are proposed as the origin of 
the transient warm dust (ISong et al..l l2005l : IWyatt et al.l 120071). One can find possible arialogo us 
transient events also in the history of our solar system (jGomes et al.ll2005l : iNesvorny et al.ll2003l ). 



In our extensive program studying debris disks around 78 F-type stars (PID: 3401, 20707, 
Moor et al., in prep., 2009) we identified four systems harboring warm disks. All four are new 
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discoveries and are located relatively close to the Sun (<60pc). Here we analyze their properties, 
discuss their possible origins and speculate about the implications on the formation and evolution 
of dust in the inner region of planetary systems. 



2. Observations and data reduction 



The four stars, HD 13246, HD 53842, HD 152598, and HD 169666 were selected for our program 
because of earlier hints for IR excess from IRAS and/or their membership in yo ung moving groups. 
HD 1 3246 and HD 53842 belong to the Tucana-Horologium moving group ([Zuckerman &: Song 
20041 ) implying an age of ~30Myr. Age indicators, such as lithium abundances and X-ray lu- 
minosities, further support this age. The galactic motion of HD 152598 (U = — 10.7, V = — 22.6, 
W = — 4.0kms~^) as well as its galactic position are consiste nt with those of th e "b3" subgroup 
of the Local Association, implying an age of 0.21ib0.07 Gyr (jAsiain et al.lll999l ). Both its frac- 
tional X-ray luminosity {ROSAT J165258. 2+314203, log{Lx/ Lboi) ~ —5.63) and the equivalent 
width of the lithium 6708 A feature (43+8 mA from spectra available in the Elodie archive. 



Moultaka et al.ll2004l'l support this age. In the case of HD 169666 isochrone fitting yields 2.1^q^ Gyr 



(jHolmberg et al.ll2007l ). Its location in the HRD indicates that this star is already leaving the main- 
sequence, supporting the conclusion of the isochrone dating. HD 169666 coincides with ROSAT 
J181909. 0+713054, its fractional X-ray luminosity is log{Lx/ Lboi) ~ —5.67. We note that this 
value is not inconsistent with the measured logfL -^/L^^;) val ues of stars with similar B — V color 



indices in the 1.9 Gyr old NGC 752 open cluster (jDaniel et al.lll994l : iGiardino et al.ll2008l ). giving 
support to the isochrone age determination. The key properties of the four stars are listed in 
Tabled! 



2.1. Spitzer observations and data reduction 



Mid- and far-IR images were obtained using the MIPS instrument (iRieke et al.ll2004l ) in small- 
field photometry mode. All four objects were measured at 24 and 70 /xm, while two of them were 
also observed at 160 fim. In addition, low-res olution spectra (R = 70-120) were obtained for each 
source using the Infrared Spectrograph (IRS. lHouck et al.ll2004l ). All these observations belong to 
our program PID 20707 except the IRS spectrum of HD 13246, which was taken from the Spitzer 
Archive (PID: 241, PI: C. Chen). In order to improve the signal-to-noise ratio and search for possible 
variability we performed follow-up IRS measurements for HD 13246 and HD 169666, obtained 41 
and 35 months after the first epochs, respectively (PID: 50538, PI: A. Kospal). 



MIPS: We started the data processing with the BCD files (pipeline version S16.1). Additional 
corrections, including flat fi eld and a background matching at 24 /xm, were performed using MOPEX 
(jMakovoz &: Marleaul l2005l ) . At 70 /xm column mean subtraction and time filtering were applied 
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following iGordon et al.l (|2007l ). BCD data were coadded with the MOPEX software. Output 
mosaics had pixels with sizes of 2'.'5, 4" and 8" at 24, 70, and 160 /im, respectively. Modified version 
of the IDLPHOT routines were used to detect sources on the final maps and to perform aperture 
photometry (including aperture correction). All of our sources were detected at 24//m, and all but 
HD 53842 at 70/xm. Both 160 /im observations provided upper limits. In the case of non-detections 
the positions derived from the 24/im images were used at 70 and 160/im for the position of the 
photometry. Photometric results are summarized in Table O 



IRS: The observations were performed in standard staring mode using high-accuracy blue peak- 
up. The data reduction began with intermediate droopres products of the SSC pipeline (version 
S15.3 for the earlier four spectra; and S18.02 and S18.1 for the follow-up obs ervations of HP 169 666 
and HD 13246, respectively). We then used the SMART reduction package (jHigdon et al.ll2004l'). in 



comb ination with IDL routines developed for the FEPS Spitzer Science Legacy program (jMever et al 



2006ll. The data red uction steps (including error budget calculation) are outlined i n detail in 



Bouwman et al. ( 20081) and have been successfully used in a series of publications (e.g. Apai et al 



2005; 



Bouwman et al 



2006; 



Pascucci et al.ll2008l . \20m . 



Results 



3.1. Spectral energy distributions and infrared excesses 



We compiled the spectral energy distribution of each star by combining IR fluxes from the IRS 
and MIPS observations with published optical and near-IR photometric data {Hipparcos, Tycho2, 
2MASS catalogs) . In Figure [1] we present the IR dat a, as well as stellar phot ospheric predictions 
obtained by fitting an ATLAS9 atmosphere model (ICastelli &: Kurucd l2003l ) to the optical and 
near-IR measurements. Metallicity data from lHolmberg et al.l (|2007l ) and computed surface gravity 
values, used as inputs to the fits, are listed in Table [TJ Since from distance est imates based on 
Hipparcos parallax data (Ivan Leeuwenll2007l ) all our stars lie in the Local Bubble (jLallement et al. 
20031 ). we assumed the extinction to be negligible. The fitted effective temperatures are included 
in Table [H Predicted photospheric fluxes for the MIPS bands are listed in Table [2l All our sources 
exhibit significant (> 3a) excess over the predicted photosphere at 24/xm. At 70/im HD 53842 was 
not detected, while the measured fluxes of the three other sources were found to be in excess. 



The IRS spectra show that the excesses are present even at shorter wavelength as they begin 
to depart from the expected photosphere at ~9, 14, 14 and 8fim in the case of HD 13246, HD 53842, 
HD 152598 and HD 169666, respectively. 



- 5 - 



3.2. Spectral features in the IRS spectra 



Figure [2] presents the IRS spectrum of HD 169666. Several broad spectral features at around 
11.3, 16.4, 23.7, 27.5, and 33.8 /im can be recognized. For comparison we plotted the spectrum of 



HD 69830 and comet Hal e- Bopp, which also exhibit features at similar wavelengths ([Beichman et al 
20051: ICrovisier et al.lll997l). These p eaks, especially the 11.3 fiui one, can be attributed to forsterite 



(jKoike et al.ll2003l : iJager et al.lll998l ). indicating the presence of micron-sized silicate grains in the 
disk of HD 169666. The spectra of HD 13246, HD 53842, and HD 152598 do not show any spectral 
features exceeding the noise. The lack of features in these disks may imply the depletion of small 
grains. 



3.3. Search for variability 

We compared the two spectra of HD 169666 obtained at different epochs. Both the continuum 
level and the spectral features were found to be remarkably unchanged. The flux differences were 
typically in the order of l%-2% over the whole spectral range. A similar comparison for HD 13246 
showed more pronounced differences especially between 10 and 20/im. However, these two measure- 
ments were processed with different pipeline versions, namely with S15.3 and S18.1. We performed 
an analysis of four calibration measurements of rj^ Dor processed with the same pipeline versions 
as HD 13246. The results indicated that the differences observed at HD 13246 could be explained 
by systematic calibrational effects related to the two pipeline versions. Thus, we conclude that no 
significant spectral change can be detected in our measurements. 



3.4. ModeUng the observed infrared excess 



We assume that the detected IR excess originates from optically thin dust confined into a 
circumstellar ring. The excess above the predicted phot osphere was fitted by a single temperature 
blackbody (as e.g. IWyatt et al.ll2005l : iRhee et al.l l2008l ). For the fitting process the IRS spectra 
were robustly averaged in five bins centered at 8, 13, 18, 24, and 32^m (the last bin was not used 
in the case of HD 169666 because of the presence of a relatively strong feature). For HD 13246 
and HD 169666, where two IRS spectra are available, we used the earlier ones because they were 
obtained closer in time to the MIPS observatio ns. An iterative m ethod was used to compute and 
apply color corrections for the MIPS data (see lMoor et al.ll2006l ). The derived dust temperatures 
{Tdust) are listed in Table [T] and the best fitting models are overplotted in Figured) 

Assuming tha t the emitting grains act l ike blackbody we estimated the ring radius using the 
following formula (iBackman &: Parescelll993l ): 



dust 

AU 



star 



dust 



(1) 
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The fractional dust luminosities of the four systems were calculated as fd = L^ust I ^hoi ^-nd are in 
the range of 0.3 — 1.8 x 10~^. The derived disk properties are listed in Table [TJ Note, that using 
the second epoch spectra for HD 13246 and HD 169666 the fitted disk parameters are not altered 
significantly. 



Discussion 



4.1. Comparison with steady state evolutionary models 



In all four systems the inferred radii of dust rings correspond to the region of the asteroid 
belt and the gas giants in our solar system. A comparison of the PR and the collision timescales 
for our dust rings (Equations (14) and (15) in iBackman fc Parescd Il993l ) implies that the grain 
evolution is dominated by collisions. The derived short lifetimes with respect to the ages of the stars 
demonstrate the second generation ('debris') n ature of the dust gr ains, and indicates the presence 
of planetesimals co-located with the dust rings. IWvatt et al.l (|2007l ) used an analytical steady state 
model to study the collisional evolution of planetesimal belts and argue that at any given age 
there is a maximum disk mass and fractional luminosity, since initially more massive disks consume 
their mass faster. We computed the maximum fractional luminosity value (fmax) for each system 



taking into account its properties as listed in Table [U and using Equation (20) in IWyatt et al 
([20071), also adopting their fixed model parameters (belt width: dr/r = 0.5, planetesimal strength: 
= 200 J kg^^, planetesimal eccentricity: e = 0.05, diameter of largest planetesimal in cascade: 
Dc = 2000 km). Comparing these values with the calculated fractional luminosities we find that 
the disk properties of HD 13246, HD 53842, and HD 152598 are consistent with a steady-state 
evolutionary scenario within the uncertainties of the model. In contrast, in the case of HD 169666 
a ratio of f /fmax > 100 implies a high dust content around this old star, inconsistently with a 
steady-state model. 



4.2. Origin of debris dust in individual systems 



HD 13246 and HD 53842 belong to the Tuc-Hor association, implying an approximate age of 30 Myr. 
In such young stars the disks are thought to be the sites of intense planet formation. In regions 
where the growing size of the largest planetesimals reach '^2000 km these protoplanets stir their 
surroundings leading to catastrophic collisions among the left-over planetesimals. A collisional 
cascade begins and produces large amounts of debris dust. Since the formation of giant plan- 
etesimals takes longer at l arger radii the site of debris production propagate outwards with time 
(|Kenvon BromlevI 12002 . booj ). 



Using the evolutionary model of iKenyon &: Bromley! (120041 ) we calculated the radius of such 
a ring for an age of 30 Myr, as a function of initial disk surface density. We found that the 3.5- 
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5AU radii derived for our disks can be explained by disk surface density of S < O.S^mmsn , 
where T^mmsn is the surface density of the minimum solar mass nebula. This parameter implies 
a relatively low initial disk mass provided that the disk is really in this evolutionary phase. After 
this peak the average dust production drops, but there may be episodic spikes due to individual 
collisions between left-over rocky protoplanets. If the dust rings around HD 13246 and HD 53842 
represent these later phases of the evolution, the initial surface density could have been higher. 

HD 152598 is older than the previous two stars, and at its age of 0.21Gyr most likely represent 
a later evolutionary phase, when dust grains are already generated by colliding minor bodies in a 
steady planetesimal ring. 



HD 169666 is a new member of a distinguished, small group of FGK-type stars harboring warm 
debris disks of unusually high fractional luminosity. Only a handful of such systenis haye been iden- 



tified so far: BD+20 307, HP 23514, HD 69830, HP 7 2905 and r/ Corvi (jSong et al. .11200,51 : lEhee et al 



20081 : iBeichman et al.ll2005l . 120061 : IWyatt et al.ll2005l ). The peculiar brightnes s of these disks is usu- 
ally explained by an episodic increase in dust production (jWvatt et al.l 120071 ). It is also a common 
property of HP 169666 and most other members of this group that they are old (>300Myr) and 
their mid-IR spectra exhibit spectral features attributed to small silicate dust grains. Such small 
grains move and evolve primarily under the influence of radiation pressure and collisions. 

In order to draw a picture on the dust evolution we computed the maximum size of dust grains 
blown out of the HP 169666 system due to the radiation pressure. Approximating the absorp- 
tion+scattering cross section by the geometric cross secti on and adopting a density of 2.5gcm~'^ 
and an albedo of 0.1, we found a 'blowout-size' of 1.6 //m (jBurns et al.lll979l ). Under the influence 
of radiation pressure alone these small grains would be removed very fast: first-order estimates 
of the grain removal from the mid-IR emitting region (T> 130 K) suggest typical timescales of 
a few months to a year. With the most prominent opacity source removed so rapidly, significant 
changes in the mid-IR spectral features are to be expected on similarly short timescales. However, 
the mid-IR spectrum of HP 169666 showed no significant variations in the course of three years. 
The constant flux levels of the spectral features indicate that the small grains are continuously 
replenished on timescales of a few years or less. The existence of warm - small or large - dust in 
this syst em covers an even long er interval. Convolving our IRS spectrum with 25 /zm filter profile 
of IRAS (jBeichman et al.lll988l ) and comparing this syn thetic fiux value of 112 mJy with the IRAS 
measurement in 1983 (101ifcl3 mJv. iMoshir et al.lll989l ). we could not see any variation within the 
measurement uncertainties. This finding suggests that the transient phase of HP 169666 lasts at 
least for a quarter of a century. It is an open question whether the small grain production is 
active also on this longer timescale. Nevertheless, the fact that most old transient systems (see 
above) exhibit spectral features might suggest that small grains are present, and their production 
mechanism is active over a significant fraction of the transient phase. Because no multi-epoch mea- 
surements exist for the other five warm disks, possible variations of the spectral features cannot be 
excluded. Petails of the variability depend on the actual mechanism of small grain replenishment 
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in the system. 



Three possible mechanisms have been put forward to explain the origin of dust in systems like 
HD 169666. First, brea kup of a large planete simal similar to the event leading to the birth of the 
Veritas asteroid family (INesvornv et al.ll2003ll. T he probability of such an event, however, is rela- 
tively low at an age of ~2 Gyr (IWyatt et al.ll2007l ). Second, dust may originate from the cont inuous 
evaporation of a super-comet as proposed for the case of HD 69830 (jBeichman et al.ll2005l ). The 
dust temperature in the disk of HD 169666 exceeds 110 K, the sublimation temperature of comets, 
consistently with this scenario. The third mechanism would be that the observed high fractional 
luminosity of these disks may be the result of an event similar to the Late Heavy Bombardment 
(|Wvatt et al.l 120071 ). Modeling the expected brightness and spectral evolution of transient disks 
around older stars in the three scenarios and comparing them with multi-epoch observations may 
help to identify the processes responsible for replenishing the observed massive dust belts in these 
extrasolar planetary systems. 
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Fig. 1. — Spectral energy distribution of the four F-type stars harboring warm debris disks. In the 
cases of HD 13246 and HD 169666 IRS spectrum of the first epochs from 2005 were plotted. 
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Fig. 2. — IRS spectrum of HD 169666 obtained on 2005 July 2. The spectrum was smoothed 
with boxcar average of 5 pixels and the photospheric contribution was subtracted. For comparison 
Spitzer spectrum of H P 69830 (jBeichman et al.l 120051 ) and the ISO spectrum of comet Hale-Bopp 
(jCrovisier et al.l 119971') are displayed. The vertical dotted lines mark the laboratory peak wave- 
lengths of forsterite (jKoike et al.ll2003l ). 
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Table 1: Stellar and disk properties 











STELLAR PROPERTIES 








DISK PROPERTIES 




Source ID 


V 


SpT 


D 


Tcff 


L. 


M. 


logg 


[Fe/H] 


Age 


Tdust 


Rdust 


f 


fmax 




(mag) 




(pc) 


(K) 


(Lo) 


(Mo) 


(cms—-') 


(dex) 


(Myr) 


[K] 


[AU] 


[10-4] 


[10-4] 


HD 13246=^ 


7.50 


F8V 


44.2 


6140 


1.55 


1.06 


4.4 


-0.17 


~30 


166±18 


3.5±0.9 


1.7 


0.8 


HD 53842 


7.46 


F5V 


56.0 


6430 


2.6 


1.20 


4.3 


-0.16 


~30 


151±24 


5.4±1.4 


0.53 


1.5 


HD 152598'= 


5.34 


FOV 


29.2 


7100 


4.8 


1.43 


4.3 


-0.21 


210±70 


135±11 


9.3±1.5 


0.35 


0.5 


HD 169666 


6.68 


F5 


53.2 


6540 


4.7 


1.35 


4.1 


-0.04 


21Q0+200 


198±13 


4.2±0.6 


1.84 


0.008 



"HD 13246 has a late type (K5) companion, CD-60 416 (HD 13246B), with a separation of ~52" on the sky. 

''HD 152598 is identified as a binary system in the CCDM catalog i Dommanget fc Nvs 20021 ). However, according to 
the NOMAD catalog, the proper motion of the proposed companion (CCDM J16530-I-3142B) deviates significantly 
from the proper motion of HD 152598, implying that a physical connection between these stars is unlikely. 



Table 2: Predicted photospheric fluxes and Spitzer/MIPS photometry 



Source ID 


P24 


F24 


P70 


F70 


P160 


F160 








[mjy] 






HD 13246 


23.4 


46.2±1.9 


2.6 


20.1±3.3 


0.5 


<60.6 


HD 53842 


21.2 


30.1±1.3 


2.3 


<27.2 


0.5 


<123.2 


HD 152598 


106.5 


134.8±5.5 


11.6 


44.3±5.0 






HD 169666 


40.2 


89.2±3.7 


4.4 


21.8±4.0 







Note. — Predicted (P24, P70, P160) and measured (F24, F70, F160) flux density values given in mJy. The 
final uncertainties were computed by adding quadratically the individual internal errors and the absolute calibration 
uncertainties. Following the MIPS Data Handbook, we adopted 4% and 7% calibration uncertainty at 24 and 70/im. 
For sources not detected at a specific wavelength, upper limits were computed as the flux measured in the source 
aperture -f 3a (all measured fluxes were positive). 



